Dynamic Nuclear Polarization in Silicon Microparticles 
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We report record high 29 Si spin polarization obtained using dynamic nuclear polarization in mi- 
crocrystalline silicon powder. Unpaired electrons in this silicon powder are due to dangling bonds in 
the amorphous region of this intrinsically heterogeneous sample. 29 Si nuclei in the amorphous region 
become polarized by forced electron-nuclear spin flips driven by off-resonant microwave radiation 
while nuclei in the crystalline region are polarized by spin diffusion across crystalline boundaries. 
Hyperpolarized silicon microparticles have long T\ relaxation times and could be used as tracers for 
magnetic resonance imaging. 

PACS numbers: 76.70.Fz, 81.07.Bc, 03.67.Lx, 82.56.-b 
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Crystalline and amorphous silicon are both of great 
technological importance. Single-crystal silicon is a work- 
horse of the electronic industry while amorphous silicon 
is widely used for photovoltaic devices. Recently, sev- 
eral promising silicon-based architectures have been pro- 
posed for quantum computing [l], 0, E|. The long co- 
herence times of 29 Si and 31 P nuclear spins in lightly 
doped silicon crystals make them attractive candidates 
for quantum bits (qubits) 0, Q • One of the major limita- 
tions of using nuclear spins as qubits is the highly mixed 
state of the nuclear spin system at thermal equilibrium 
S @1- Several methods can be used to hyperpolar- 
ize nuclear spin systems. For instance, optical pumping 
techniques have been very successful in studies of elec- 
tron spin physics in gallium arsenide quantum wells and 
quantum dots @, [l(| but so far have produced compara- 
tively low nuclear spin polarization in silicon [HI . [l3| . 
Microwave induced dynamic nuclear polarization (DNP) 
is a powerful technique, where an enhancement of sev- 
eral orders of magnitude of the nuclear spin polarization 
can be generated by microwave irradiation of paramag- 
netic centers coupled to the nuclear spin system through 
the hyperfine interaction 1J, |15| . Since the discovery of 
DNP in 1953 lri 17 1, there have been several reports of 
3 but all of them were done 



DNP in silicon 1 



either at high temperatures or low magnetic fields lim- 
iting the 29 Si spin polarization to a few percent It is 



imperative for the emerging fields of spintronics 21| and 
nanotechnology to develop methods to achieve high 29 Si 
spin polarization in both bulk and nanoscale samples. 

Here, we report the first DNP experiments in undoped 
silicon microparticles. The amorphous region of these 
particles has a very high concentration of structural de- 



-300- 



-400 



-500 



— DNP 4hours 
-■ DNP 100s 
-- No DNP 4 hour wait 




-2 -4 -6 
Frequency (kHz) 



-10 



'Author to whom correspondence should be addressed. Electronic 
address : anatolydCSphysics. harvard. edu 



FIG. 1: 29 Si NMR spectra acquired at 1.4K for different mi- 
crowave irradiation times in B=2.35T (f =19.89MHz). The 
spectrum acquired without the microwave irradiation is shown 
using the dashed line. 



fccts - dangling bonds [22|]. 29 Si nuclei in this region are 
polarized by microwave irradiation of the hyperfine cou- 
pled electron -nuclear spin system. The density of para- 
magnetic impurities in the crystalline cores is low which 
results in the long nuclear relaxation times. Nuclei in 
those regions are polarized by nuclear spin diffusion [23[ 
across crystalline boundaries. Several orders of magni- 
tude enhancement of the NMR signal combined with the 
long Ti relaxation time makes this a valuable technique 
for hyperpolarizing silicon nanoparticles, with potential 
applications in medical imaging QUI and NMR studies 
of nanoscale systems (e.g. silicon nanowires). 

Our sample is a polycrystalline powder obtained from 
Alfa Aesar. It is 99.999% pure with particle sizes dis- 
tributed between 1 and 5 /im. According to X-ray diffrac- 
tion analysis of this sample 80% of the sample is amor- 
phous while 20% is crystalline and individual crystal- 
lites are larger than 200nm. Our DNP experiments were 
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performed at 1.4 K in a 2.35 T superconducting NMR 
magnet. An Oxford NMR Spectrostat was operated 
in single-shot mode to cool the sample down to 1.4 K. 
The microwave source was a 90 mW Gunn diode source 
(Millitcch). The NMR spectrometer used was a Bruker 
Avance system with a home-built probe, containing a 
horn antenna for the microwave irradiation of the sam- 
ple, and a solenoidal NMR rf coil, that were in direct 
contact with pumped liquid helium. Details of a similar 
probe design have been described elsewhere (26j|. 

Figure 1 shows the 29 Si NMR spectra acquired without 
and with microwave irradiation at 66.25 GHz for different 
irradiation times. For short irradiation times the NMR 
spectrum is inhomogeneously broadened with a 2 kHz 
linewidth (full width at half of a maximum height) due 
to the distribution of distances between nuclear and elec- 
tron spins in the amorphous region. The sharp feature in 
the spectrum is due to nuclei in the crystalline cores and 
appears for microwave irradiation times longer than an 
hour. The 300 Hz linewidth of this feature is consistent 
with previous NMR measurements on silicon crystalline 
powder samples 27] and can be explained by a combi- 
nation of dipole-dipole interaction between 29 Si nuclei 
arranged in the crystalline lattice and magnetic suscep- 
tibility broadening. The average DNP enhancement of 
the nuclear polarization for the amorphous phase inferred 
from Fig.l is 150, corresponding to a 29 Si polarization of 
5%. We believe this to be the highest 29 Si polarization 
yet achieved. Nuclear polarization of the crystalline re- 
gion after 4 h of DNP is found to be ~ 2%. This estimate 
comes from the ratio of the integrated NMR signals from 
amorphous and crystalline regions and takes into account 
that crystalline region is only 20% of the sample. 

We performed ESR studies to characterize the nature 
of the paramagnetic centers in our sample. The g-factor 
value g — 2.006, obtained from our ESR spectra is con- 
sistent with earlier measurements in amorphous silicon 
[22] ] that identified the defects as dangling bonds. Figure 
2(a) shows the low temperature ESR spectrum which was 
acquired in the rapid-passage regime with the spectrom- 
eter tuned to dispersive response. The electron spin Ti 
relaxation time was estimated to be approximately 30 /is 
using the rapid-passage phase method (28l |. 

Figure 2(b) shows the dependence of the NMR signal 
amplitude on the microwave frequency for a microwave 
irradiation time of 200 s acquired at 2.35T and 1.4K. 
It is characteristic of the DNP thermal mixing mecha- 
nism. This is a typical DNP mechanism in insulators 
with a high concentration of paramagnetic impurities at 
low temperatures [3, [ID, [2^] . Microwave irradiation just 
below the ESR central frequency cools the electron spin 
dipolar reservoir while microwave irradiation just above 
the ESR central frequency creates a negative tempera- 
ture for the electron spin dipolar reservoir. Since the 
ESR linewidth is larger than the nuclear Larmor fre- 
quency, the electron spin flip-flop interaction mediated 
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FIG. 2: (a) The ESR spectrum acquired at 3.6 K and 9.27 
GHz. (b) The dependence of the NMR signal amplitude on 
microwave frequency acquired at 2.35T and 1.4K for a mi- 
crowave irradiation time of 200 s. 



by nuclear spin flips equalizes the electron spin dipolar 
temperature with the nuclear spin temperature. 

The microwave irradiation time dependence of the 
NMR signal amplitude is shown in Fig. 3(a). There are 
two distinct time scales. Within the first hour 29 Si nuclei 
in the amorphous region are hyperpolarized through the 
thermal mixing mechanism. Beyond the first hour the 
NMR signal continues to increase because of spin diffu- 
sion into the crystalline cores. There is still an increase of 
the signal for times after 4 h. This is consistent with very 
slow spin diffusion in the network of naturally abundant 
29 Si spins mediated by XY terms of the dipolar interac- 
tion: 
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ratio for 29 Si). The vector between spins i and j, fy, 
satisfies • z — cos 9ij . Magnetization transport is 
described by the diffusion equation: 



dM 



= DAM 



(2) 



where D 



-^jr is the diffusion constant [23], T 2 is the 
nuclear transverse relaxation time and b is the average 
distance between neighboring 29 Si nuclei. Assuming in 
our case T2 ~ 5.6ms 27 1 and b « 0.5571 -1 / 3 w 0.41nm 



12 3 4 

Microwave irradiation time (hours) 





0- 






Z3 




.a 


-100- 


3. 




as 


-200- 






CO 




rr 


-300- 


5 




z 




CO 


-400- 



4000 -4000 -8000 
Frequency (Hz) 

FIG. 3: (a) The dependence of the NMR signal amplitude 
on microwave irradiation time. The solid line is to guide the 
eye. (b) Two 29 Si NMR spectra, the solid line is the spectrum 
acquired after 4 h of DNP at 1.4K with a 5° tipping pulse and 
the dashed line is the spectrum acquired after the microwave 
irradiation was turned off and the sample was warmed up to 
240K (a 90° tipping pulse was used). 



[30J, where n is the density of 29 Si nuclei, we estimate 
D ^ 6 x 10" 15 cm 2 /s. 

Once the nuclear polarization builds up inside the crys- 
talline cores of the particles it can be stored for a time 
on the order of the nuclear Ti relaxation time. For ex- 
ample, Fig 3(b) shows two 29 Si NMR spectra, one (solid 
line) was acquired after 4 h of DNP at 1.4K with a short 
(5°) NMR pulse, while the other one (dashed line) was ac- 
quired after the microwave irradiation was turned off and 
the sample was slowly warmed up to 240K. The broad 
line from the amorphous phase quickly disappeared due 
to the fast relaxation rate of the amorphous domain and 
only the narrow component from the crystalline cores 
survived after an hour long warming up procedure. The 
resulted 29 Si polarization is about 3000 times higher than 
the equilibrium polarization at 240K. 

In order to estimate the 29 Si nuclear Ti in the crys- 
talline part of the sample, we prepared nuclear spin po- 
larization by DNP at 1.4K for 5 h, allowing for nuclear 
spin diffusion to penetrate inside the crystalline cores. 
Afterwards, the microwave radiation was turned off and 
the cryostat was warmed up to 4.2K to enable the con- 
tinuous flow operation. Figure 4(a) shows the decay of 
the NMR signal amplitude from the crystalline cores to- 
gether with a simulation (solid line) using the diffusion 
equation (Eq.|2J) for a spherical particle with a radius of 
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FIG. 4: (a) The decay of the NMR signal from crystalline sil- 
icon after DNP for 5 h at 1.4K. We used 15° NMR pulses to 
monitor the signal decay. The solid line is the calculation of 
the decay due to the spin diffusion of nuclear polarization to 
the surface, (b) Simulated nuclear spin magnetization den- 
sity profiles in a spherical particle with a radius of 350nm 
corresponding to different time points in (a): right after DNP 
(labeled hours), I hour and 5 hours after DNP. 



350nm. In combination with the large initial polariza- 
tion, it is the long Ti observed here that underlies the 
potential use of these particles as tracers. 

In the simulation we set the magnetization density on 
the surface of the particle to a finite value for the 5 h DNP 
period and to a zero right afterwards. The nuclear spin 
magnetization density profiles obtained using this simula- 
tion are shown in Fig. 4(b). The agreement between the 
experiment and the simulation is satisfactory given the 
simplicity of the model. The decay of the spin magneti- 
zation in this model depends on the particle size and the 
nuclear spin diffusion coefficient. Although the decay is 
nonexponential, we calculate the nuclear spin-lattice re- 
laxation time Ti as the time when magnetization decays 
to 1/e of the initial value. It scales as Ti oc where 
R is the radius of a particle. There is an interplay be- 
tween the initial magnetization and the Ti on one side, 
and the isotopic abundance of 29 Si on the other. High 
29 Si concentration leads to a larger initial magnetization 
but also to a larger nuclear spin diffusion coefficient and 
thus a shorter Ti. The result is that isotopic enrichment 
reduces the storage time. Figure 5 shows simulations 
of the 29 Si magnetization decay due to spin diffusion and 
instantaneous relaxation at the surface for different parti- 
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curity Agency (NSA) under Army Research Office (ARO) 
contract number DAAD190310125, the NSF and DARPA 
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FIG. 5: (a) Simulations of the 29 Si magnetization decay due 
to spin diffusion and instantaneous relaxation at the surface 
for particles with 300nm radius, (b) Simulated Ti relaxation 
time as a function of the particle size. Solid line is for the nat- 
urally abundant 29 Si and dashed line is for the 29 Si enriched 
sample. 



cle sizes and for naturally abundant (4.67%) and enriched 
(100%) 29 Si. The ratio of spin diffusion coefficients for 
enriched and naturally abundant 29 Si is about 10 31 1. 
Although the initial magnetization density is 21.4 times 
larger for enriched silicon, the magnetization decay is also 
much faster due to faster diffusion and there is always a 
crossing between the two curves at long times. 

In summary, we have extended previous DNP exper- 
iments in silicon to lower temperature and higher mag- 
netic field and have achieved record high 29 Si spin po- 
larization in silicon microparticles. This technique might 
also be used for polarizing nuclei in thin films of crys- 
talline silicon if a layer of amorphous silicon is created 
on top (e.g., by ion implantation). Nuclear spin diffusion 



in Si abundant crystals is significantly faster which will 
make it possible to polarize thicker films. The decay of 
nuclear spin polarization in the crystalline part of the 
sample is due to nuclear spin diffusion to the surface and 
subsequent relaxation by paramagnetic impurities. The 
use of hyperpolarized nanoparticles as tracers depends on 
both the size and time scales of interest in a given appli- 
cation. While moderately sized particles can retain their 
polarization for a relatively long time, the relaxation be- 
comes very fast for the smallest nanoparticles. 
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